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ABSTRACT 



We present here the analysis of about 19,500 new star hours of low ecliptic 
latitude observations (\b\ < 20°) obtained by the Hubble Space Telescope's Fine 
Guidance Sensors over a time span of more than n ine years; which is an ad dition 
to the ~ 12, 000 star hours previously analyzed by ISchlichting et al.l (120091 ). Our 
search for stellar occultations by small Kuiper belt objects (KBOs) yielded one 
new candidate event corresponding to a body with a 530 ± 70m radius at a 
distance of about 40 AU. Using bootstrap simulations, we estimate a probability 
of ~ 5%, that this event is due to random statistical fluctuations within the new 
dat a set. Combining this n ew event with the single KBO occultation reported 
by ISchlichting et al.l ( 120091 ) we arrive at the following results: 1) The ecliptic 
latitudes of 6.6° and 14.4° of the two events are consistent with the observed 
inclination distribution of larger, lOOkm-sized KBOs. 2) Assuming that small, 
sub-km sized KBOs have the same ecliptic latitude distribution as their larger 
counterparts, we find an ecliptic surface density of KBOs with radii larger than 
250 m of N(r > 250 m) = x 10 7 deg~ 2 ; if sub-km sized KBOs have instead 

a uniform ecliptic latitude distribution for —20° < b < 20° then N(r > 250 m) = 
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4.4+2'8 x 10 6 deg~ 2 . This is the best measurement of the surface density of sub-km 
sized KBOs to date. 3) Assuming the KBO size distribution can be well described 
by a single power law given by iV(> r) oc r 1 ~ 9 , where N(> r) is the number of 
KBOs with radii greater than r, and q is the power law index, we find q = 3.8±0.2 
and q = 3.6±0.2 for a KBO ecliptic latitude distribution that follows the observed 
distribution for larger, 100-km sized KBOs and a uniform KBO ecliptic latitude 
distribution for —20° < b < 20°, respectively. 4) Regardless of the exact power 
law, our results suggest that small KBOs are numerous enough to satisfy the 
required supply rate for the Jupiter family comets. 5) We can rule out a single 
power law below the break with q > 4.0 at 2a, confirming a strong deficit of 
sub-km sized KBOs compared to a population extrapolated from objects with 
r > 45 km. This suggests that small KBOs are undergoing collisional erosion 
and that the Kuiper belt is a true analogue to the dust producing debris disks 
observed around other stars. 

Subject headings: Occultations — Kuiper belt: general — Comets: general - 
Planets and satellites: formation — Methods: observational — Techniques: pho- 
tometric 



1. INTRODUCTION 

The Kuiper belt consists of a disk of icy objects and is located at the outskirts of our 
planetary system just beyon d the orbit of Nep tune. Since the discovery of the first Kuiper 
belt object (KBO) in 1992 (IJewitt et al.lll992l ) well over 1000 objects have been detected. 
Dedic ated surveys hay e been revealing an intricate dynamical structure of the Kuiper belt 



(e.g., iPetit et al.ll201ll ). have been investigating the multiplicity of KBOs (e.g.. iNo 



20081 ). and have been measuring the Kuiper belt size-distribution (e.g.. lBernstein et al 



et al. 



200J). 



In the Kuiper belt, planet formation never reached completion because runaway growth of 
the planetary embryos was interrupted, most likely, by an increase in velocity dispersion 
of the planetesimals. The size distribution of larger KBOs provides a snap-shot of planet 
formation that was erased elsewhere in the solar system, where planet formation went all 
the way to completion. 

The cumulative size distribution of KBOs larger than about 45 km in radius is well 
described by a single power-law given by 



N(> r) oc r 1_ 



T 



where N(> r) is the number of objects with radii greater than r, and q is the power-law index. 
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Observations find that the power-law index for large KBOs is q ~ 4.5 (IFuentes &: Holman 
20081 ; iFraser et al.ll2008l ). Since this size distribution is a relic of the accretion history in the 
Kuiper belt, it provides val uable insights into the formation of large KBOs a nd the planet 
formation process itself (e.g. ISternlll996l ; iDavis fc Farinellalll997l ; lKenyonll2002l ). Coagulation 
models of planet formation reproduce the observed size distr ibution of large KBOs well (e.g., 
Kenvon fc Luulll999[ lKenvonll2002[ ISchlichting fc Sarilboilh . 



Observatio ns suggest that there is a break in the power-law size distribution at smaller 
KBO sizes (e.g. Bemstein et al]|2004l ; iFuentes fc Holman|[2008l ; IFraser et aDbood ; ISchlichting et al 



2009: Fuentes et al 



201~oL The break in the size distribution is generally attributed to col- 



lisions that br eak-up small KBOs (i.e.,r < 45 km) and s u bsequently modif y their size dis- 



tribution (e.g. iDohnanyil Il969t iKenyon fc Bromley! 12004 ; |Pan &: Saril 120051 ). If this is so 



then the location of the break (i.e., the break radius) constrains the time period over which 
destructive collisions have been occurring in the Kuiper belt. This is because with time, 
collisions move the break radius to ever larger sizes. Furthermore, the power-law index be- 
low the break radius constrains the material properties of KBOs. Small KBOs that are in 
collisional equilibrium and that are held together predominantly by material strength are 
expected to have a size distribution with power-law index q = 3.5, which is the so-called 
Dohnanyi spectrum (IDohnanyil Il969l ). The size distribution will be shallower with a power- 
law index of q ~ 3 if small KBOs are held t ogether predominantly by gravity (i.e., they 
are effectively rubble piles) (IPan fc Saril 120051 ). The break radius and the power-law index 
below the break, therefore, constrain the collisional history of the Kuiper belt and reveal the 
material properties of small KBOs, respectively. 

Our interest in studying the abundance and size distribution of small (r < 10km) KBOs 
is threefold. First of all, we want to determine their overall abundance to test if they are 
numerous enough to supply the Jupiter Family comets. Second, we want to measure their 
size distribution in order to constrain collisional evolution in the Kuiper belt and establish it 
as a true analog of the dust producing debris disk observed around other stars. And third, 
we would like to constrain the small KBO size distribution well enough to infer their material 
properties as discussed above. 

Unfortunately, not much is known about the KBO size distribution significantly below 
the break, because objects with radii less than ~ 10km are too faint to be detected in 
reflected ligh t . These objects can, however, b e detected indirectly by stellar occultations 
(jBaileyi Il976l ; lDysonlll992t lAxelrod et al.lll992l ). A small KBO crossing the line of sight of 
a star will partially obscure the stellar light, whi ch, under suitable circumstances, can be 



detected in the star's light curve. Both opt i cal fe.g. iRooues et al 



Bianco et al.l 120091 ; ISchlichting et al.l 120091 ; iBianco et al 



2010 



2006 



Bickerton et al. 2008 



Wang et al.l 120101 ) and X-ray 
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( jChang et al.ll2006U201lh occultation s earches have been conducted to probe the population 
in the Kuiper belt. IChang et al.l ( 120061 ) searched for occultation signatures in arc hival RXTE 



X-ray data of Scorpius-Xl and reported a surprisingly high event rate. However, I Jones et al. 



( 120081 ) showed that most of the dips in the Scorpius-Xl light curves are artificial effects 
caused by the RXTE photomultiplier . Only 12 of the original 58 events were not ruled out 
as artifacts; however iBickerton et al.l ( 120081 ) point out that the duration and/or depth of 
most of the 12 remain i ng events are inconsiste n t with the diffraction signature of occultation 



events. 



Rogues et al.l (120061 ) , IBickerton et al.l ( 120081 ) , and IWang et al.l (120101 ) all conducted 
occultation surveys in the optical regime. None of these surveys reported any detections 
of objects in the Kuiper belt, which is not surprising given the low expected event rate of 
these surveys. Ground-based observations may suffer from a high rate of false positives due 
to atmospheric scintillation. The only ground-based system that is currently able to screen 
events caused by atmospheric scintillation and other kinds of interference is the Taiwanese 
American Occultation Survey (TA OS), which consists of four telescopes that observe the 
same position simultaneously (e.g.. lAlcock et al.l 120031 ) . TAOS collected over 500,000 star 



hours with 4 or 5 Hz sampling frequen c y and reported no detections so far (IBianco et al. 



2010l ). Recently iDoressoundiram et al.l ( 120121 ) analyzed about 140,000 star hours sampled 
at 1 HZ by the COROT satellite and reported up to 15 possible KBO occultation events. 
The first measurement of the surface density of sub-km sized KBOs has been reported by 
Schlichting et al.l (|2009|). who analyzed 12,00 star hours of archival HST data with 40 Hz 
sampling frequency. ISchlichting et al.l (120091 ) reported one detection and found an ecliptic 
surface density of KBOs with radii larger than 250 m of 2.lti'y x 10 7 deg -2 . 



In this paper we present the analysis of new data obtained by the Fine Guidance Sensors 
(FGSs) on board the Hubble Space Telescope (HST) over a time span of more than nine 
years. We searched this data set for stellar occultations by small KBOs with the aim to 
measure their abundance, to constrain their size distribution and to address the question of 
whether they are ubiquitous enough to satisfy the supply needed to explain the observed 
frequency of Jupiter family comets. 

This paper is structured as follows. We describe our HST-FGS survey in section [2j De- 
tails about our data analysis and detection efficiency follow in sections [3] and IU respectively. 
In section [5] we present the results from our survey. We calculate the abundance of sub-km 
sized KBOs, compare our findings with other Kuiper belt surveys and conclude in section |6j 
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THE HST-FGS OCCULTATION SURVEY 



The HST Fine Guidance Sensors (FGSs) are an integral part of the HST pointing control 
system, ensuring a pointing stability of HST at the milli-arc-second level over exposure times 
of tens of minutes. Each FGS is a dual-axis white light shearing interferometer and the 
photon count of each FGS is recoded by 4 PMTs (e.g. Fine Guidance Sensor Instrument 



Handbook: http://www.stsci.edu/hst/fgs/documents/instrumenthandbook/fgs_ihb.pdf). In 
addition to ensuring the pointing stability of HST, the FGSs have been used as science 
instruments. For example, the FGS was used as an astrometric in strument to measure 
the precise distances to eight galactic Cepheids ([Benedict et al. as an interferometer 

to resolve binary O-stars in the Carina Nebula ( jNelan et al.l 120041 ). and as a high speed 
photometer to study transit light c urves of the planet orbiting HD and asteroseismology of 
the star itself rtGilliland et allboilh . 



For over a decade, the three FGSs have been collecting a large number of photometric 
measurements of the guide stars with 40 Hz time resolution. In this paper we search more 
than nine years of this archival data set for photometric signatures consistent with the 
occultation of the guide star by a small object in the Kuiper belt. The 40 Hz sampling 
frequency, the large number of star hours of FGS observations that have been obtained over 
the lifetime of HST, and the stability inherent to space based observations make this data 
set ideal for the search of serendipitous stellar occultations by small objects in the outer 
solar system. 

The en tire data set, including l ow and high ecliptic latitude observations and the data 
analyzed in ISchlichting et al.l (120091 ) . consists of ~ 90,000 star hours. The signal-to-noise 
ratio, defined here as the square root of the mean number of photon counts of all 4 PMTs 
combined in a 40 Hz interval, as a function of star hours is shown in Figure [1] and ranges 
from about 10 to 100, reflecting the 9 < V < 14 range in magnitude of the HST guide stars. 
Nominal HST operation uses two FGSs for guiding, with each FGS observing its own guide 
star. We are therefore able to remove false positives due to instrumental effects and, for 
example, day/night-time variations due to HST's orbit around the Earth. Figure [2] shows 
the HST integration time as a function of ecliptic latitude. About 35% of the observations 
are taken in the low ecliptic latitude region (i.e., —20 deg < i < +20 deg) of the Kuiper Belt. 
The remaining high-ecliptic latitude observations provide an excellent control sample which 
we use to test for possible false positives. Figure [3] shows the distribution of angular sizes 
of the HST guide stars in our data set in units of the Fresnel angular scale at 40 AU and for 
observations at a wavelength of 600 nm. The Fresnel angular scale is given by ~ a/A/2o, 
where A is the wavelength of the observations and a the semi-major axis of the Kuiper belt. 
About 63 % of the stars in our data set subtend angular sizes less than 0.5 Fresnel angular 
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Fig. 1. — Distribution of star hours as a function of the mean signal-to-noise ratio, S/N, in 
a 40 Hz interval for the entire FGS data set. The signal-to- noise ratio is defined here as the 
square root of the mean number of photon counts in a 40 Hz interval. 
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Fig. 2. — Distribution of star hours as a function of ecliptic latitude for the FGS data set. 



scales at a distance of 40 AU, and 84 % of all stars in our data set have angular sizes of less 
than one Fresnel scale. The diffraction pattern that is produced by a sub-km sized KBO 
occulting an extended background star is smoothed over the finite stellar disk. This effect 
becomes clearly noticeable for stars that subtend sizes larger than about 0.5 Fresnel scales 
and reduces the detectability of occultation events around such stars. The distribution of 
the finite angular sizes of the stars are taken into account when we calculate the detection 
efficiency of our survey in section HI 

Finally we analyzed the entire data set for correlated noise by calculating the auto- 
correlation function with lags of 0.025 seconds. Most of the data sets, each consisting of 
the observations obtained over one HST orbit, are free of statistical significant correlated 
noise. About 17% of the roughly 120,000 the data sets contained correlated noise, which is 
significant at, or above, the 4.5<r level for a given data set and significant at, or above, the la 
level for the analysis of the entire data, which consisted of 120,000 independent trials. This 
correlated noise is often due to slopes, i.e. long-term variability, in the data. Such long-term 
trends may affect the results from a bootstrap analysis, however, we show that the two data 
files containing candidate events for which we use a bootstrap analysis to estimate their 
significance do not have any statistical significant correlated noise (see section 5 for details). 



-8 - 




Fig. 3. — Distribution of angular sizes of the HST guide stars in our data set, in units of the 
angular Fresnel scale. The angular sizes of the guide stars were calculated by fitting 2MASS 
JHK and USN0-B1 BR photometry with a black-body spectrum. 
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3. DATA ANALYSIS 

The 40 Hz time resolution of our survey allows for the detection of the actual diffrac- 
tion pattern rather than a simple decrease in the photon counts. Our detection algorithm 
therefore employs a template search using theoretical light curves and performs \ 2 fitting of 
the templates to the data. This template fitting procedure improves the sensitivity of the 
survey compared to algorithms that search only for dips in the light curve and aids with the 
identification of false positives. 

In general, the occultation diffraction pattern is determined by the size of the KBO, the 
angular size of the star, the wavelength range of the observations and the impact parameter 
between the star and the KBO. The typical KBO sizes that our survey is most likely to detect 
is determined by the signal-to- noise of our data, because small KBOs are more numerous 
than their larger counterparts and the amplitude of the occultation signal scales as the cross- 
section of the KBO. Given a typical signal-to-noise ratio of 15 in a 40 Hz interval (see Figure 
[IJ, our survey is most likely to detect KBO occultation events caused by objects that are a 
few hundred meters in radius. Since these objects are significantly smaller than the Fresnel 
scale, they give rise to occultation events that are in the Fraunhofer diffraction regime. This 
implies that the diffraction pattern itself is not sensitive to the exact shape of the KBO, but 
the amplitude of the diffraction pattern scales linearly as the cross section of the occulter. 
This significantly reduces the number of templates that need to be implemented in the search 
algorithm. Furthermore, the templates that we use in our search algorithm treat the stars 
as point sources. The reason for this is twofold: First, a large fraction (~ 63%) of the stars 
in our survey have angular sizes below 0.5 Fresnel scales, which implies that the point source 
template is a good approximation for these stars. Second, using templates with finite angular 
size stars hurts rather than helps the detection efficiency because the typical uncertainties 
in the estimated angular radii are too large (i.e. they are typically about 40%). Since 
the diffraction pattern is wavelength dependent, we integrate the light curve templates that 
we use in our detection algorithm over a wavelength range of the FGS observations which 
extends from 400 to 700 nm. Finally, our search algorithm includes light curve templates 
for impact parameters between the KBO and background star ranging from to 1.0 Fresnel 
scales in steps of 0.2 Fresnel scales. 

For a given impact parameter, our theoretical light curves have three free parameters. 
The first is the mean number of photon counts, m, which is the normalization of the light 
curve. The second is the amplitude of the occultation, A, which is proportional to the cross 
section of the KBO, and the third is the duration of the occultation, which is inversely 
proportional to the relative velocities between HST and the KBO. We obtain values for the 
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first two parameters from our data by minimizing \ 2 where 



2 _ (di - modelj) 



and 

model; = m(A(L - 1) + 1) (3) 

where Zj is our theoretical light curve template. We can simply solve for the values of A and m 
from our data and therefore only need to perform a full search the third free parameter, which 
is the duration of the occultation and time. The duration of the occultation is independent 
of the object size, and mainly determined by the ratio of the Fresnel scale to the relative 
speed between HST and the KBO perpendicular to the line of sight. This relative speed is 
given by the combination of HST's velocity around the Earth, Earth's velocity around the 
Sun and the velocity of the KBO itself. We use this information to restrict the parameter 
space for the template widths in our search such that we are sensitive to KBOs located at the 
distance of the Kuiper belt between 30 AU and 60 AU, but allow for KBO random velocities 

Of Up to fRepler = 5km/s. 

Our search algorithm identifies occultation candidates by calculating their A% 2 , which 
we define as the difference between the x 2 value calculated for a horizontal line, corresponding 
to no event, and the x 2 value derived from best fit light curve template. Occultation events 
have large A% 2 since they are poorly fitted by a constant straight line, but well matched by 
the light curve template. If the noise properties were identical over the entire data set, then 
the probability that a given occultation candidate is due to random noise can be characterized 
by a single value of A% 2 for all observations. In reality however, the noise properties are 
different from observation to observation; especially non-Poisson tails in the photon counts 
distribution will give rise to slightly different Ax 2 distributions. Therefore, ideally, we would 
determine a unique detection criterion for each individual HST orbit. However, this would 
require to simulate each data set, which contains about an hour of observations in a single 
HST orbit, over the entire length of our survey, which is not feasible due to the enormous 
computational resources that would be required. Instead, we perform bootstrap simulations 
over all the FGS data sets toge ther and us e this to estimate the typical A% 2 value that we 



use in our detection algorithm ( lEfronlll982l ). Using bootstrap simulations for estimating the 
significance of candidate events is justified as long as there is no correlated noise in the data. 
For all occultation candidates that exceed this detection threshold, we determined their 
statistical significance, i.e. the probability that they are due to random noise, by extensive 
bootstrap simulations of the individual data sets. Our detection algorithm flagged all events 
for which the template fit of the diffraction pattern was better than 15a and that had a 
Ax 2 > 63. This detection criterion corresponds to about 0.5 false-positive detections over 
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the 19,500 new star hours of low ecliptic latitude observations that are analyzed in this paper. 
All flagged events that were solely due to one single low or high point were ignored in the 
data analysis, but included as false-positives when calculating the significance of candidate 
occultation events. This eliminates a large number of otherwise flagged events that may, for 
example, be caused by cosmic rays. 



4. DETECTION EFFICENCY 

The ability to detect an occultation event of a given size KBO depends on the impact 
parameter of the KBO, the duration of the event, the angular size of the star and the signal- 
to-noise ratio of the data. We calculate the detection efficiency of our survey by planting 
synthetic events with different radii into the FGS data and analyzing this modified data set 
with the same search algorithm that we used to analyze the original FGS data with the 
same significance threshold of Ax 2 > 63. The synthetic events correspond to KBO sizes 
ranging from 200 m < r < 850 m, they have impact parameters from to 2.5 Fresnel scales 
and a relative velocity distribution that is identical to that of the actual FGS observations. 
To account for the finite angular sizes of the stars we generated light curve templates with 
stellar angular radii of 0.1, 0.3, 0.4, 0.6, 0.8 and 1 Fresnel scales. The detection efficiency 
of our survey was calculated using the angular size distribution of the FGS guide stars as 
shown in Figure [3] and the synthetic events were implanted in a subset of the FGS data that 
had the same signal-to-noise properties as the whole data set shown in Figure [TJ Figure 0] 
shows the detection efficiency, i](r), as a function of KBO radius. We normalize our detection 
efficiency for a given size KBO to 1 for an effective detection cross section with a radius of 
one Fresnel scale. The detection efficiency of our survey is ~ 0.04 for objects with r = 200 m 
and ~ 0.75 for KBOs with r = 400 m located at 40 AU. 



5. RESULTS 

5.1. One New Occultation Candidate 

Among the ~ 40 candidate events that were flagged with a Ax 2 > 63 in the new data 
set, all but one turned out to be false-positives (see subsection 15.21 for a detailed discussion 
of the false-positives), leaving us with one new occultation candidate event. Figure [5] shows 
the candidate event with the best fit template from our search algorithm. The red points 
and and error bars represent the FGS data with Poisson error bars. The actual noise in 
this observation is about 8% larger than Poisson noise, which is due to additional noise 
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Fig. 4. — Detection efficency, ^(r), as a function of KBO radius. The detection efficiency is 
normalized by an impact parameter equal to 1 Fresnel scale, i.e., a detection efficiency of 1 
means that the effective cross-section for detection has a radius equal to one Fresnel scale. 
The detection efficiency of our survey is ~ 4% for objects with r = 200 m and 75% for KBOs 
with r = 400 m. 
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sources such as dark counts, which contribute about 3 to 6 counts for a given PMT in a 40Hz 
interval. The mean signal-to-noise ratio in a 40 Hz interval for this HST orbit of observations 
is ~ 10. The best fit x 2 /dof from our detection algorithm is 27.3/28. Each FGS provides 
two independent PMT readings and we confirmed that the occultation signature is present 
in both of the these independent photon counts. The position of the star is R.A. =64. 74065°, 
DEC=28.13064° (J2000), which translates to an ecliptic latitude of +6.6°. We obtained a 
spectrum of this star on 1 2 February 2012 wit h the Echelle Spectrograph and Imager (ESI) 



on the Keck II telescope (ISheinis et al.ll2002l ). We used ESI's echellete mode with a 0.5" 



slit width, providing spectral coverage from 3900 to 11000 Angstroms simultaneously with a 
resolving power (A/AA) of 8100. We exposed for 300 seconds, achieving a median signal-to- 
noise of 75 between 6000 and 9 000 Angstroms. Analy sis of the stellar spectrum and fitting 



the JHK bands from 2MASS (ISkrutskie et al.ll2006l ) yields an stellar angular radius and 



effective temperature of 0.58 ± 0.06 Fresnel scales and ~ 5000 K, respectively. Using our 
best estimate for the stellar angular radius, we find that the best fit parameters yield a 
KBO size of r = 530 ± 70m and, assuming a circular orbit, a distance of 35 ± 9AU. We 
note here, that within the uncertainties of the actual stellar angular radius, that light curve 
templates with smaller angular radii give a somewhat better fit to the data than templates 
with larger stellar radii. For objects an circular orbits around the sun, two solutions can fit 
the duration of the event. However, the second solution corresponds to a distance of 0.2 AU 
from the Earth and an objects size of ~ 50 m, and is therefore unlikely. It is also unlikely 
that the occulting object was located in the Asteroid belt, since the expected occultation 
rate from Asteroids is about two orders of magnitude less than our implied rate. In addition, 
an Asteroid would have to have an eccentricity of 0.2 or greater to match the duration of 
the observed occultation candidate. 

We estimate the probability that this candidate event is due to statistical fluctuations 
using bootstrap simulations. This approach is justified as long as there is no correlated 
noise in the data. We calculate the autocorrelation for lags between and 1 seconds for 
the HST orbit of observations that contained the event and find that the autocorrelation 
function is consistent with zero (see Figure [7]). Using the data from the HST orbit, that 
contained the event, we removed the event itself and simulated 3.1 x 10 6 star hours, which 
corresponds to 161 times the low ecliptic latitude observations analyzed in this paper. This 
calculation required ~ 2300 CPU days of computing power. Figure [8] shows the cumulative 
number of false-positives, Nf_ p as a function of A% 2 . The number of false-positives was 
normalized to 19,500 star hours, which correspond to the length of the entire low ecliptic 
latitude observations analyzed here. In the entire bootstrap analysis we obtained 8 events 
with a Ax 2 > 71.9. This implies a probability of ~ 5% that events like the occultation 
candidate with A% 2 = 71.9 are caused by random statistical fluctuations over the entire low 
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Fig. 5. — Photon counts as function of time of the candidate occultation event observed 
by FGS3. The red points and error bars are the FGS data points with Poisson error bars, 
the dashed blue line is the theoretical light curve, and the blue squares correspond to the 
theoretical light curve template in our detection algorithm integrated over 40Hz intervals. 
We note here that the actual noise for this observation is about 8% larger than poison noise, 
due to additional noise sources such as, for example, dark counts, which contribute about 3 
to 6 counts per PMT in a 40Hz Interval. The best fit x 2 /dof from our detection algorithm 
is 27.3/28. The star has an ecliptic latitude of +6.6° and its angular radius and effective 
temperature are 0.58 ± 0.06 Fresnel scales and ~ 5000 K, respectively. The position of the 
star is R.A.=64.74065°, DEC=28. 13064° (J2000) and its estimated V-magnitude is 13.9. 
Assuming a circular orbit, the best fit parameters yield a KBO size of r = 530 ± 70m and a 
distance of 35 ± 9AU. 
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Fig. 6. — Photon counts as a function of time during the occultation event observed by FGS3. 
The red-solid line represents the 40 Hz photon counts from FGS3, while the black-dashed 
line shows the displacement of the guide star along the x-axis, Ax, in the FGS relative to its 
position at time = -2 seconds. The occultation event occurred (at time=0) approximately 
one second after HST began a commanded small angle maneuver associated with a planned 
displacement of the active science instrument's aperture on the sky. During this maneuver 
FGS3 actively tracked its guide star, keeping it at interferometric null, just as it did prior 
to the maneuver. Thus the decrease of photon counts near time = is not correlated with 
or caused by the telescope repointing. Note, the displacement along the FGS y-axis was far 
smaller and is therefore not plotted here. 
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ecliptic latitude observations analyzed in this paper. 
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Fig. 7. — Autocorrelation function for lags between and 1 seconds for the HST orbit of 
observations containing the candidate event shown in Figure [51 We note that the autocorre- 
lation at lag is, by definition, equal to one. The red-dotted lines shows the upper and lower 
la errors of the autocorrelation. The autocorrelation function of this data is consistent with 
zero. 



5.2. False-Positives 

Among the ~ 40 candidate events that were flagged because they had a A% 2 > 63, 
all but one turned out to be false-positives. The most common false-positives were due to 
what looks like a slower read out of the photon counts (see Figure [9]) or showed a correlation 
between the signature in the photon counts and the displacement of the guide star from its 
null position on the FGS (see Figure [TO]) . Figure ITOl shows an example of a false-positive that 
shows a strong correlation between the number of photon counts and the displacement of the 
guide star from its null position on the FGS. The jitter introduced due to the displacement 
of the guide star from its null position causes an up to 3% change in the photon counts and 
is therefore only detected by our search algorithm for stars that have photon counts above a 
few thousand in a single 40Hz interval. The jitter due to the guide star's displacement has 
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Fig. 8. — Cumulative number of false-positives, Nf- P as a function of A% 2 for the candidate 
event normalized to 19,500 star hours, which correspond to the length of the entire low 
ecliptic latitude observations analyzed in this paper. The false-positives were obtained from 
bootstrap simulations using the one orbit of HST observations (~ 50 minutes) in which we 
found the candidate event. From these bootstrap simulations we find a probability of ~ 5% 
that events like the occultation candidate with A% 2 = 71.9 are caused by random statistical 
fluctuations in the entire low ecliptic latitude data set. 
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Fig. 9. — Photon counts as a function of time. An example of a false-positive found in the 
FGS data, which is due to what looks like a slower read out of the photon counts. False- 
positives like this one are due to update problems that are encountered as the telemetry from 
HST goes through various transfer stations and ground stations on its way to the archive. 



- 19 - 



a characteristic frequency of 1Hz and is therefore easy to identify 
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Fig. 10. — An example of a false-positive that is due to jitter introduced by the displacement 
of the guide star from its null position. The red-solid line represents the photon counts of 
the FGS data and the black-dashed line the corresponding displacement of the guide star 
relative to its null position along the y-axis, as recorded by the FGS, as a function of time. 



5.3. High-Ecliptic latitude Control Sample 

In addition to the 19,500 star hours of low ecliptic latitude observations, |6| < 20°, we 
also analyzed ~ 36, 000 star hours of high ecliptic latitude, |6| > 20°, observations and used 
these as a control sample. We analyzed the high ecliptic latitude observations with exactly 
the same detection algorithm as the low ecliptic latitude data. A total of ~ 70 candidate 
events were flagged with Ax 2 > 63. All but one of these events were due to either jitter, 
induced due to the displacement of the guide star from its null position (see Figure [TO]) . 
or a lower read out frequency of the PMTs (see Figure [9]). The only event which was not 
caused by either of these two effects is shown in Figure [LTJ The ecliptic latitude of this 
event is 81.5° and it has a A% 2 = 72.3. Calculation of the autocorrelation function of this 
orbit of HST observations showed that there is no statistical significant correlated noise in 
the data. Using bootstrap simulations we find a 21% probability that this events is due to 
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random statistical fluctuations (see Figure IT2l . given the entire high ecliptic latitude data 
set analyzed in this paper. This implies a ~ 79% chance that this flagged event is due to 
a high inclination KBO, but this interpretation seems unlikely given the known inclination 
distribution of large KBOs. 




-0.2 -0.1 0.0 0.1 0.2 
time relative to mid eclipse[s] 

Fig. 11. — Photon counts as function of time of the only flagged event in our control sample 
that could not be contributed to either jitter, induced due to the displacement of the guide 
star from its null position (see Figure [TO]) , or a lower read out frequency of the PMTs 
(see Figure [9]). The red points and error bars are the FGS data points with Poisson error 
bars, the dashed blue line is the theoretical light curve, and the blue squares correspond to 
the theoretical light curve integrated over 40Hz intervals. The best fit x 2 /dof is 13.5/17. 
Bootstrap simulations yield a probability of ~ 21% that events like this with A% 2 = 72.3 
are caused by random statistical fluctuations over the high ecliptic latitude control sample 
analyzed in this paper. This implies a ~ 79% chance that this flagged event is due to a 
high inclination KBO, but this interpretation seems unlikely given the known inclination 
distribution of large KBOs. 
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Fig. 12. — Cumulative number of false-positives, Nf_ p as a function of A% 2 for the flagged 
event found in the high ecliptic latitude control sample. The number of false-positives is 
normalized to 36,000 star hours, which correspond to the length of the entire high ecliptic 
latitude observations analyzed in this paper. The false-positives plotted here were obtained 
from bootstrap simulations using the one orbit of HST observations in which this event 
occurred. From these bootstrap simulations we find a probability of ~ 21% that events like 
this with Ax 2 = 72.3 are caused by random statistical fluctuations over the high ecliptic 
latitude control sample analyzed in this paper. 
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Discussion & Conclusions 



We combine the one candidate KBO occultation event presented in this paper, which 
we found in 19,500 star ho urs of low ecliptic latitu d e HST-FGS observation s, with the single 
event that was reported by lSchlichting et al.l (120091 ) . ISchlichting et al.l ( 120091 ) analyzed 12,000 
star hours of low ecliptic latitude HST-FGS observations and reported one KBO occultation 
event at an ecliptic latitude of +14.4°. First we test whether the ecliptic latitudes of the two 
events are consistent with the observed inclination distribution of larger KBOs and then we 
use the two events to estimate the abundance of sub-km sized KBOs. 



6.1. Inclination Distribution 

In this subsection, we test if the observed ecliptic latitude distribution of the two oc- 
cultation events is consistent with the inclination distribution of larger KBOs inferred from 
direct searches. 

We first calculate for each star in our data set the ecliptic latitude, and the amount of 
time it was observed by the FGSs. Figure [2] shows the distribution of star hours as a function 
of ecliptic latitude. We denote this function by £fgs(/3)- Using the inclination distribution 
of KBOs, P( ^kbo), from lEUiot et al.l (120051 ). we randomly draw inclinations from P(«kbo) 
and for each random declination we choose a random ecliptic latitude sin(/3) = sin(z) sin(A) 
0. Here the ecliptic longitude, A, is a random number drawn from a uniform distribution 
between and 2tt. This give us the probability distribution of the instant aneous eclip t ic lat i- 
tude distribution, Pkbo(/3), associated with inclination distribution from lEUiot et al.l ( 120051 ) . 
Figure [T3l shows the product of Pkbo(/3) and £fgs(/3) (denoted by -Pfgs(/3)), and gives the 
probability that our survey will detect a KB O at a given ecliptic latitude, assuming the 
KBOs are drawn from the lEUiot et al.l (120051 ) inclination distribution. We assumed here 
that KBOs follow circular orbits and that the S/N distribution of the FGS data is the same 
at all ecliptic latitudes. 

In order to see if our observations are consistent with the inclination distribution of 
large KBOs, we calculate the probability distribution for observing two events at different 
latitudes using the detection probability as a function of ecliptic latitude shown in Figure [131 
The log-probability of our ecliptic latitude observations is: Y^j m -fFGs(/3j), where /3j are our 
observed ecliptic latitudes of +6.6° and +14.4°. In order to put this into context we need 
to estimate the log-probability distribution given our data set. This is done using Monte- 



1 Specifically, we use the numerical inclination distribution presented in Fig 20. of lElliot et al.l ([20051 ) 
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Carlo simulations, which consist of drawing 10,000 randoms pairs of ecliptic latitudes and 
calculating the log-probability for each pair. Figure [14] shows the probability distribution, 
where the arrow denotes the value of the log-probability of our observations. This plot 
suggests that our observed ecliptic latitudes are consistent with ecliptic latitudes drawn 
from the inclination distribution of large KBOs. However, this conclusion would change if 
the flagged event in our control sample (discussed in section [573]) was due to an actual KBO. 
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Fig. 13. — Occultaion detection probability per degree, Pfgsi as a function of ecliptic latitude 
for the low ecliptic latitude observations \b\ < 20° of the FGS data set. The detection 
probability was calculated from the ecliptic latitude distribution of the FGS data shown in 
Figure [2] Note, we assumed that the KBO ecliptic latitude distribution is symmetric about 
the ecliptic and ignored the small 1.6° inclination of the Kuiper belt plane relative to the 
ecliptic. 



6.2. Abundance of Sub-km sized KBOs 



We now combine the results from this paper with the work of 
to estimate the abundance of sub-km sized KBOs. 



Schlichting et al 



(120091 ) 
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Fig. 14. — Probability distribution derived from Monte-Carlo simulations, which consisted 
of drawing 10,000 randoms pairs of of ecliptic latitudes from the distribution shown in Figure 
[13] and calculating the probability for each pair. The arrow denotes the log-probability value 
of the of our observations, suggesting that our observed ecliptic latitudes are consistent with 
ecliptic latitudes drawn from the inclination distribution of large KBOs. 
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The number of occultation events is given by 

N events = -2v rel F I ^ I v (r)^- dN ^ 6) dbdr (4) 

•J T-min J —b 

where v re i = 22km/s is the typical relative velocity between the KBO and the observer, b 
is the ecliptic latitude, At/Ab is the time observed per degree in ecliptic latitude, as shown 
in Figure [2, and F = 1.3 km is the Fresnel scale. The sky density of KBOs is both a 
function of ecliptic latitude, b, and the KBO radius, r. Therefore, in order to estimate the 
total number of KBOs of a given size or their corresponding sky density, we need to make 
an assumption regarding their ecliptic latitude distribution. Unfortunately very little is 
currently know about the ecliptic latitude distribution of sub-km sized objects in the Kuiper 
belt. We therefore estimate the total number of KBOs for two very different ecliptic latitude 
distributions, currently both are consistent with t he ecliptic latitudes of th e candidate event 



presented in this paper and the event reported by ISchlichting et al.l (120091 ) . In the first case 
we assume that the KBO latitude distribution, /(&), is uniform between —20° and +20°, 
such that f(b) = 1 for —20° < b < 20° and zero otherwise. In the second case we assume 
that small sub-km sized KBOs follow the same ecliptic latitude distribution as their larger 



100 km -sized counterparts and use the ecliptic latitude distribution provided in Elliot et al. 



(120051 ) in Figure 14. We further assume that the KBO size distribution follows a power law, 
such that it can be written as N(r,b) = nor~ q+1 f(b) where no is the normalization factor 
for the cumulative surface density of KBOs. Substituting for dN(r, b)/dr in Equation 0] and 
solving for n we get 

^events / r \ 

n o — t. : • (5) 

2v rel F(q - 1) ! r Z: v(r)r-« dr J b _ b f(b)f b db 

Evaluating Equation [5] assuming a uniform KBO ecliptic latitude distribution for —20° < 
b < 20° yields a cumulative KBO surface density of 

N(r > 250 m) = 4.4^ j x 10 6 deg" 2 . (6) 

Similarly, evaluating Equation [51 assuming that th e small, sub-km sized KBOs follow the 



ecliptic latitude distribution from Elliot et al.l ( 120051 ). yields a cumulative KBO surface den- 



sity averaged over the ecliptic < 5°) of 

N(r > 250 m) = 1.1 j;^ x 10 7 deg" 2 . (7) 

When evaluating the integral over r in Equation [5], we assumed q — 4. We note however that 
the value for the cumulative KBO surface density at r=250 m only weakly depends on the 
exact choice for q. We therefore find an ecliptic KBO abundance for bodies with r > 250 m 
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that ranges, depending on the actual inclination distribution of sub-km sized KBOs, between 
4.4 x 10 6 deg -2 to 1.1 x 10 7 deg~ 2 . This is the best measurement of the surface density of 
sub-km sized KBOs to date and about a factor of 2 lower than the first results published by 



Schlichting et al.l (120091 ). Figure [TS] displays the results from the FGS survey and summarizes 
published upper limits from various works. The red point plotted at r=250m with the upper 
and lower error bars in Figure [151 gives the best estimate of the KBO surface density around 
the ecliptic (—5° < b < 5°) from our survey with la errors, assuming that sub-km sized 
KBOs follow the same inclination distribution as their larger, 100-km sized, cousins. The 
upper and lower red curves correspond to our upper and lower 95% confidence level, which 
are derived without assuming any size distribution. This limit and the red point would be 
a factor of 2.4 lower if sub-km sized KBOs would have an ecliptic latitude distribution that 
is close to uniform for —20° < b < 20°. The jump between 500 m and 600 m in the upper 
limit curve is due to the fact that below 500 m it is ca lculated for 2 even ts, whereas above 
600 m for no events. The 95% upper limit from TAOS faianco et alJboid ) is about a factor 



of 2 lower than the one derived from our FGS survey, if sub-km sized KBOs have the same 
ecliptic latitude distribution as their larger, lOOkm-sized, counterparts, and about the same 
if they follow a uniform ecliptic latitude distribution between —20° and +20°. 

Assuming that the KBO size distribution can be well described by a single power law that 
is parameterized by iV(> r) oc r 1 ^ q , where iV(> r) is the number of KBOs with radii greater 
than r, and q is the power law index, we can use the above estimated KBO abundances 
to calculate the power law index below the observed break in the KBO size distribution. 
Assuming a break radius of 45k m and a corresponding cumulative KBO surface density 



of 5.4deg (IFuentes et al.l 120091 ) we find q = 3.6 ± 0.2 and q = 3.8 ±0.2 for a uniform 



KBO ecliptic latitude distribution and a KBO latitude distribution from lEUiot et al.l ( 120051 ). 
respectively. 



Trilling et al.l ( 120101 ) found that, in contrast to larger KBOs with sizes above the break 
radius, fainter (R > 26) KBOs, which have sizes below the break radius, are dominated 
by dynamically excited objects (i > 5°). If this result applies all the way to sub-km sized 
KBOs, then this suggests that the inclination distribution for sub-km sized KBOs may be 
dominated by dynamically excited objects and that the true abundance of sub-km sized 
KBOs, and their corresponding size distribution power-law index, may in fact lie somewhere 
between N(r > 250 m) = 4.4 x 10 6 — 1.1 x 10 7 deg~ 2 , and q = 3.6 — 3.8, that we estimated 
above (see Figure IT5l). In addition, we can rule out a single power law below the break with 
q > 4.0 at 2a. 
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Fig. 15. — Cumulative KBO size distribution as a function of KBO radius. The results from 
our FGS survey are presented in three different ways: (i) the red point plotted at r=250 m 
with the upper and lower error bars gives the best estimate of the KBO surface density 
around the ecliptic (—5° < b < 5°) from our survey with la errors, assuming that sub-km 
sized KBOs follow the same inclination distribution as their larger, 100-km sized, cousins. 

(ii) The upper and lower red curves correspond to our upper and lower 95% confidence level. 

(iii) The la range for our best estimate of the power law size distribution index, q = 3.8 ±0.2 
i s given by the yello w shaded region normalized to N(> r) = 5.4 deg -2 at a radius of 45km 



(jFuentes et al.l 120091 1 . The plotted red point, the upper and lower limits and the power law 
size distributions were all derived assuming that sub-km sized KBOs have the same ecliptic 
latitude distribution as their larger, lOOkm-sized, counterparts. All would be a factor of 
2.4 lower if sub-km sized KBOs would have an ecliptic latitude distribution that is close to 
uniform for —20° < b < 20°. The area between the two blue horizontal lines defines the 



required scattered disk KBO abundances from IVolk fc Malhotral (120081 ) in order to supply 
the Jupiter Family comets. In addition, 95% upper limits from various optical KB O occul- 
tation surveys are sho wn as bla ck symbols with ar rows (circles IBianco et al.l (120091). s quare 
Bickerton et all ( 120081 ) . triangle IWang et al.l (120101). and diamond iRoques et al.l (120061 )) and 
as a dashed back line for TAOS (IBianco et al.ll2010l ). 



-28 - 



6.3. Conclusion 



We presented here the analysis of ~ 19, 500 new star hours of low ecliptic latitude 
< 20°) archival data that was obtained by the HST over a time span of more than 
nine years. Our search for stellar occultations by small Kuiper belt objects (KBOs) yielded 
one new candidate event, which, assuming a circular orbit, corresponds to a body with 
a ~ 500 m radius located at a distance of about 40 AU. Using bootstrap simulations, we 
estimate a probability of ~ 5%, that this event is due to random statistical fluctuations 
within the analyzed data set. Com bining this new event with the single KBO occultation 
reported by ISchlichting et al.l (120091 ) we show that their ecliptic latitudes of 6.6° and 14.4°, 
respectively, are consistent with the observed inclination distribution of larger, 100-km sized 
KBOs. 



Assuming that the new candidate event and the event previously reported by ISchlichting et al 



(120091 ) are indeed genuine KBO occultations and that small, sub-km sized KBOs have the 
same ecliptic latitude distribution as larger KBOs, we find an ecliptic surface density of 
KBOs with radii larger than 250m of 1.1^0 7 x 10 7 deg -2 . If sub-km sized KBOs have in- 
stead a uniform ecliptic latitude distribution for —20° < b < 20° we find N(r > 250 m) = 
4.4I2 8 x 10 6 deg -2 . The ecliptic latitudes of the two events are consistent with both a uniform 
ecl iptic latitud e distr ibution and the ecliptic latitude distribution of larger KBOs published 
by lEiliot et al.l (120051 ) . These estimated KBO abundances provide the best measurements of 
the abundance of sub-km sized KBOs to date and are, although consistent wit hin la, about 
a factor of 2 lower than previous results published by ISchlichting et al.l (120091 ) . 



Assuming that the KBO size distribution for bodies with radii smaller than the break 
radius can be well described by a single power law that is parameterized by iV(> r) oc r 1-9 , 
where N(> r) is the number of KBOs with radii greater than r, and q is the power law index, 
we find q = 3.6 ± 0.2 and q = 3.8 ± 0.2 for a uniform KBO ecliptic latitude distribution and 
a KBO ecliptic latitude distribution that follows the observed distribution for large KBOs, 
respectively. These results are consistent with a power-law index of q = 3.5, correspond ing to 
a collisional cascade consisting of material strength dominated bodies (jDohnanyilll969l ) that 
all have the same constant velocity dispersion, within better than 1 and 2cr, respectively. 
We caution however that the actual size distribution of small KBOs is likely to exhibit 
significant deviations from a single power law due to possible changes from gravity to martial 
stren gth dominated bodies and waves that may e xist in the small KBO size distribution 



(e.g. lO'Brien fc Greenbergl 120031 ; iPan fc Sarill2005l ). In addition, the velocity dispersion of 
the bodies in the cascade may not be constant and instead evolve as a function of size, 
which results in a size distributions that is significantly steeper than the one derived without 
velocity evolution (e.g. the standard q = 3.5 power-law index of the Dohnanyi differential 
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size spectrum can change to an index as large as q = 4) (IPan Schlichtingi 120121 ) . 

Finally, our findings suggest that small KBOs are numero us enough to sati s fy the re- 
quired supply rate for the Jupiter family comets calculated by IVolk &: Malhotral (120081 ) for 
scattered disk objects. We can rule out a single power law below the break with q > 4.0 at 
2cr confirming a strong deficit of sub-km sized KBOs compared to a population extrapolated 
from objects with r > 45 km. This suggests that small KBOs are undergoing collisional ero- 
sion and that the Kuiper belt is a true analogue to the dust producing debris disks observed 
around other stars. 
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